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SUMMARY 

In an earlier paper, we discussed the probability that the determinant of a matrix undergoes the least 
change upon perturbation of one of its elements, provided that most or all of the elements of the 
matrix are chosen at random and that the randomly chosen elements have a fixed probability of being 
non-zero. In this paper, we derive approximate formulas for that probability by assuming that the 
terms in the permanent of a matrix are independent of one another, and we apply that assumption 
to several classes of matrices. In the course of deriving those formulas, we identified several integer 
sequences that are not listed on Sloane's Web site. Preprint submitted to Numer. Linear Algebra 
Appl. on June 28, 2007 
Prepared using niaauth.cis 
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1. Introduction 

In an earlier paper [1] , we discussed the problem of finding the probabihty that the determinant 
of a matrix undergoes the least change under perturbation of one of its elements. In this paper, 
we consider only the case where the randomly chosen matrix elements are values of a continuous 
(real) random variable, and we derive approximate formulas for that probability by assuming 
that the terms in the permanent of a matrix are mutually independent. Denote the determinant 
of any matrix A by detA. For an (n + 1) x (n + 1) matrix M„+i, the expansion of detM„+i 
via row i is 

det M„+i = rriiiMii H h m^jMij H h mj„+iM„+i, (1) 

where is the element of Mn+i at the intersection of row i and column j, and My is the 
cofactor of rriij. Mij can be written as (— I)*"'"-' det S",!, where Sn is the n x n submatrix of 
Mn+i which is obtained by deleting row i and column j. Thus the probability that detM„_|_i 
undergoes the least change upon perturbation of element rriy is equal to the probability that 
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I det Sn I is as small as possible. As in the earlier paper, we treat the following three classes of 

n X n matrices Sn- 

(i) Matrices An in which all the elements are values of mutually independent random 
variables each of which has probability r of being non-zero and probability 1 — r of 
being 0, where < r < 1. 

(ii) Matrices i?„ in which all but one of the diagonal elements are set to 1 (i.e., bu = 1 
for i 1, where bu is the special diagonal element), and 6ii and all the off-diagonal 
elements are as in (i). Thus 
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(iii) Matrices C„ in which all the diagonal elements are set to 1 and all the off-diagonal 
elements are as in (i). Thus 
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The matrices 5„ described in (i), (ii), and (iii) above will be said to be matrices of type A„, B„, 
and C„, respectively. Every randomly chosen element of a matrix of any of these types (i.e., 
every element of a matrix of type A„, every off-diagonal element of a matrix of type B„ or C„, 
and the special diagonal element 611 of a matrix of type B„) will be called a variable element; 
the remaining elements will be called fixeA elements. If X is a continuous set such that G X 
(as is assumed here), and X is a random variable whose values are selected from X, then the 
probability function P for X is defined by 



P{X = 0) 
P{X ^ 0) 



1 -r 
r 



Since X is continuous, P can be expressed in terms of a density function q as 

DCX^ P{X eD)= I q{x)dx, 
Jd 

where 

^ q{x) ^1, / q{x)dx = 1 — r, / q{x)dx = r 
J{0} Jx-{0} 

If q is continuous on X — {0} (as is assumed here), then 

1 — r, x = 



(3) 



(4) 



P{X = x) = 



0, 



a; G X - {0} 



(5) 
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As found in [1], the probability that dctM„+i undergoes the least change upon perturbation 
of element rriij is equal to the probability that det Sn = u, where m = if iS^ is of type A„ or 
B„, and m = 1 if 5„ is of type C„. In [1], the probability that det5„ = u was formulated in 
terms of binary matrices (matrices of O's and I's), as follows: For a matrix S'„ of type A„, B„, 
or C„, Sn is the n x n binary matrix with Sij = 1 if ^ 0, and = if Sjj = 0. Moreover, 
Sij is a variable (resp. fixed) element of Sn if Sij is a variable (resp. fixed) element of Sn, and 
Sn is of type A„ (resp. B„, C„) if Sn is of type A„ (resp. B„, C„). Every variable element 
of Sn has probability r of being 1, and probability 1 — r of being 0. Denote the permanent of 
any matrix A by per A. The expansion of the permanent of Sn via the permutation group 6„ 
on the set {1, 2, 3, . . . , n} is 

perS'„= ^ S<7(l)lSa(2)2---S<7(n)n (6) 

Since Sn is a binary matrix, every term in this expansion is either or 1, so the value of per Sn 
is a non-negative integer. In [1], it was shown that the probability that det Sn = u is equal 
to the probability that per Sn = u, and that Sn has non-zero (positive) probability of having 
permanent u if and only if every term in the expansion of per 5„ via (3„ that contains at least 
one variable element is 0. Thus if Sn is of type A„ or B„, then Sn has non-zero probability of 
having determinant if and only if every term in the expansion of per Sn is 0; and if Sn is of 
type C„, then Sn has non-zero probability of having determinant 1 if and only if the only non- 
zero term in the expansion of per Sn is the "diagonal" term, C11C22C33 • • • c„„. A matrix of type 
An can be converted to a matrix of type B„ by replacing all but one of the diagonal elements 
with a fixed element 1. The latter matrix can in turn be converted to a matrix of type C„ by 
replacing the sole variable element on the main diagonal (element 611) with a fixed element 
1. Therefore, we might expect to find that there exist analytic formulas for -PperB„=o(^) 
terms of PperA„=o(^)' and Pperc„=i(^) in ^rms of PperB„=o(»')> where Ppers„=«(^) denotes 
the probability that perS^„ = u. In [1], the exact probability -PperS„=u('') ^ specific type of 
matrix (A„, B„, or C„) was formulated in terms of the numbers of matrices 5„ of that type 
which have i variable elements with a value of 1 in the expansion of their permanent (where i 
ranged from to some Zmax that depended on the type of matrix). In this paper, we will derive 
an approximate probability Qp^^g that per Sn = u, by assuming that all the terms in 

the expansion of per Sn via 6„ are independent of one another. For example, the expansion 
of the permanent of a matrix of type A3 is 

perA3 = aiia22C33 -I- 013032021 -f- 012023031 + 013031022 
+011023032 + 012621033 

If 023 = 0, the terms 012023031 and 011623032 must both be zero, but we will ignore that 
relationship. Instead, for every m < n and each of the three types of matrices (A„, B„, and 
C„), we will consider only the number of terms in the expansion of the permanent of a matrix 
of that type which have m variable elements, regardless of the values of those variable elements 
or the connections between different terms in the expansion. 
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2. Formulation of (3p„s„=„(r-) 



There are n! terms Sct(i)iSo-(2)2 • ' • Sa{n)n in (6), because the number of elements a G S„ is n!. 
For each term, let m be the mimber of variable elements it contains, and let En{m) be the 
number of terms with m variable elements. Then 

n 

m— 1 

where r™ is the probability that all the variable elements in a term with m variable elements 
are non-zero, 1 — r"* is the probability that at least one variable element in a term with m 
variable elements is 0, and (1 — 7'™)-E„(m) jg ^-^^^ probability that all the terms with m variable 
elements are 0. What remains to be done is to determine £'„(to) for every m < n and each of 
the three types of matrices. 

2.1. Type An 

For type A„, every term in the expansion of per A„ has n variable elements, so En{m) = for 
every m < n. Since there are n! terms, we obtain 



OperA„=oW = (l-r-T 



(7) 



2.2. TypeCn 

For type C„, let W„(m) denote En{m), so that 



Proposition 2.1. 



m=l 

Ty„(n-l) = n-W„_i(n-l) 



(8) 

(9) 



Proof Every term in the expansion of per C„ that has n — 1 variable elements contains just 
one fixed (diagonal) element (namely, ca for some i with 1 <i<n). For every i, there exists 
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Figure 1. Cross shapes deleted from Cn to create the (n — 1) x (n — 1) submatrices of C„ 

a natural one-to-one correspondence between the set Tu of terms in the expansion of per C„ 
that contain element cu (and have n — 1 variable elements of C„ apiece) and the set of 
terms in the expansion of per (5^ that have n — 1 variable elements each, where (5^ is the 
(n — 1) X (n — 1) submatrix of (?„ which is formed by deleting the "cross shape" that consists 
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of row i and column i. There arc n such cross shapes (one for each diagonal element of C^), 

as shown in Fig. 1; moreover, for every i, 

tildeCl^ is a matrix of type C„_i. Thus we obtain (9). ■ 



Proposition 2.2. For every m with 1 ^ m ^ n — 1, 



Wn{m) = ^ • Wmim) (10) 
mi 



Proof Let m, k be such that 1 < m < n — 1 and k = n — m. A term in the expansion 
of per Cn has m variable elements if and only if it has k fixed (diagonal) elements (namely, 
Ci^i^,. . .,Ci^i^: for some ii,...,ik with 1 < ii < i2 < ■ ■ ■ < ik < n). For every fc-tuple {h,.. .,ik) 
with 1 < ii < Z2 < • • • < ife < n, there exists a natural one-to-one correspondence between the 
set Ti^^i^^ of terms in the expansion of perC„ that contain elements Cjiij, . . . ,Ci^i^ (and 
have m variable elements of Cn apiece) and the set 12 {j^^ j^j. of terms in the expansion of 

per C'ri*^' "'**"^ that have m variable elements each, where Cn^''"'^''^ is the m x m submatrix 
of Cn which is formed by deleting the cross shapes that correspond to ii,...,ik (i.e., the 
submatrix of C„ which is formed by deleting rows ii,...,ik and columns ii,...,ik)- The 
number of fc-tuples (ii, . . . , ik) with 1 < ii < ^2 < ■ • • < ifc < n is „Cfe; moreover, for each such 
/c-tuple, 

tildeCn^''"'^''^ is a matrix of type Cm- Also, „Cfc = " ™ , since k = n — m and „C 

m! 



n^n—m 



{n — m)\m\ ml 



Thus we obtain (10). 



Proposition 2.3. For n> 1, 



n-l 



^^n(n) = ^(-l)"-^+i-„P,-i (11) 



Proof Wn(n) is the number of permutations of the numbers 1, 2, 3, . . . , n with no fixed points 
(i.e., the number of permutations that have no cycles of length one). Such permutations are 
called derangements. In 1713, Pierre de Montmort [2] proved that 



W„(n)=n-W„_i(n-l) + (-l)" (12) 
5 
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Thus we obtain (11) recursively as follows: 



Wn{n) = n- + (13) 

= „Pi-P^„_i(n-l) + (-l)"-„Po 

= n • [(n - l)Wn-2{n - 2) + (-1)""'] + (-1)" 

= n{n - 1) • W„_2(n - 2) + n • (-1)""^ + (-1)" 

2 

= „P2 • Wn-2{n - 2) + ^(-l)"-^+i . „P^_i 

1=1 



fe 



= • • • = „Pfe ■ W„_fe(n - fc) + ^(-1)""^ + ' • nPj-l 
= • • • = nPn-1 ■ W,{1) + ^(-1)"-J + 1 . „P^_i 

^■ — 1 



n-1 

= ^(-l)"-^ + l-„P,_i 

i=i 

The last step follows from the fact that M^i(l) = 0. 
Proposition 2.4. For every m with < m < n, 



f-lV 



Proof The only term in the expansion of per C„ that has no variable elements is Cn • • • c„ 
(-1)0 

so W„(0) = 1 = „Po ■ This proves that (14) holds for m = 0. By (11), 



Wn{n) = ^(-l)"-^"+i-„P,_i 

n— 1 I 

-i+i . 



,.=1 (n-i + 1)! 
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This proves that (14) holds for m = n (the last step follows from the fact that — Tj — 
1 — 1 = 0), so we can assume that 1 < m < n — 1. By (10) and (11), 
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The values of W„(m) for n = 1,...,6 (and m = 0, ...,n) are given in Table I. The 







Table I. Values of Wn{m) for n 
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sequence {W„(n)}„>o is number A000166 in Sloane's list [3]. The first (n = 0) term, 
which has no meaning for our purposes, is given there as 1. For n > 1, Wn{n) is not only 
the number of derangements of the numbers 1,2,3, ...,n, but also the permanent of the 
n X n binary matrix with O's on the main diagonal and I's everywhere else. The sequences 
{W„(2)}„>2, {W„(3)}„>3, {T4^„(4)}„>4, and {W„(5)}„>5 arc given in Sloane's list (numbers 
A00217, A007290, A060008, and A060836, respectively). There is no m > 5 for which the 
sequence {Wn{m)}n>m is given in that list. An alternative method of determining the value 
of W„{m,) is to use the cycle structure of the permutation group (3„. For n = 5, the types 
and numbers of cycles associated with elements of 65 that are applicable to each value of m 
are listed in Table II, together with the values of W5(m) computed from them. For m = 4 
and m = 5, there are two different types of cycles each. The permutations that correspond to 
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m = 4 (which index the terms in the expansion of per Sn that have 4 variable elements) are 
those that can be expressed as a 4-cycle (such as (1234)) and those that can be expressed as 
a product of two 2-cycles (such as (12)(34)). Similarly, the permutations that correspond to 
m = 5 (which index the terms with 5 variable elements) are the one that can be expressed 
as the 5-cycle (12345) and those that can be expressed as a product of one 2-cycle and one 
3-cycle (such as (12)(345)). 



Table II. Computation of W5{m) from group tabic for ©5 
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2.3. TypeBn 

For type B„, let Vn{m) denote En{m), so that 
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Proposition 2.5. For every m with 1 < m < n, 
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Proof Let 611 be the sole variable element on the main diagonal of a matrix of type B„. If 
we delete the L-shape that consists of row 1 and column 1, we obtain a matrix of type C„_i. 
Therefore, 

Vn{m) - Wn-i{m - 1) = Wn{m) - Wn-i{m) (17) 



By (14), 
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The values of Vn{m) for n = 1, . . . , 8 (and m = 0, . . . ,n) are given in Table III. The sequence 
{Ki}n>o is number A000255 in Sloane's list [3]. The first (n = 0) term, which has no meaning 
for our purposes, is given there as 1. For n > 1, V^(n) is the number of permutations of the 
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Table III. Values of Vn{m) for n = 1, . . . , 8 (and m = 0, . . . ,n) 
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numbers 1,2,3,. ..,n+l such that there is no element k which is mapped to fc+ 1; it is also the 
permanent of the n x n binary matrix that has O's in all but one of the elements on the main 

diagonal, and I's everywhere else. The only value of m for which the sequence {Vn{m)}n>m is 
given in Sloane's list is 3 (where the sequence {V„(3)}„>3 is number A045943). 



3. Comparison of Formulas 
For n = 3, the formulas for the approximate probabilities <3perS3=«('') 

Qpe.A3=oW = {i-rr- = {i-rr 



m=l 

3 

2 



m=l 

The formulas for the exact probabihties P^^j. §^^^{r) derived in [1] are 

^peri3=oW = (l-r)9 + 9r(l-r)«+36r2(l-rf + 78r3(l-r)6 

+90r''(l - rf + 45r'5(l - r)^ + 6r^{l - rf 
^perS3=oW = (l-r)^ + 6r(l-rf + 13r2(l-r)5 + 10r3(l-r)4 

^perC3=iW = (l-rf +6r{l-rf + 12r^{l-r)^ + 6r^{l-rf 

All six of the functions given above (for n = 3) are graphed vs. r in Fig. 2. Also, Fig. 3 shows 
the six curves for n = 5 vs. r. For every n > 1 and all three types of matrices (A„, 
tilde^n, and C„), 

n 

Y,En{m) = n\ (18) 

m=0 

Thus the differences in the values of En{m) for the different types of matrices stem entirely 
from differences in the distribution of n\ over the individual values of m. For type A„, that 
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Figure 2. Graphs of QperS3=«('') and P^^.g^^^i'^) ^s- r 
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Figure 3. Graphs of Q^^.s^^ui'^) and Ppe,s,=„(r-) vs. r 

distribution is trivial (namely, En{n) ~ n\, and En{m) — for m < n). In the case of type 
C„, we found an analytic formula for the distribution: 

M^„(m) =„P„, m==0,l,2,...,n 

As illustrated earlier, the distribution for type C„ can also be derived from the cycle structure 
of the permutation group 6„. Moreover, we have shown that the distribution for type B„ can be 
derived from that for type C„, by using the relation Vn{m)—Wn-i{m—l) = Wn{m)—Wn-i{m) . 
As established in [1], there is an analytic formula for the coefficients in P^^^q =i(^)i since 
they can be derived from characteristics of acyclic digraphs with vertex set {1, 2,3,..., n}. 
For several values of n, the coefficients in ^'per a„=o('') ^^'^ -PpcrB„=o('') were also presented 
in that paper, but they were determined by computer, using an algorithm that took a 
binary matrix of a given type as input and computed its permanent. It remains to be 
seen whether there are explicit formulas for the coefficients in ^'per a„=o('') -fperB„=o(^)' 
and whether there is an analytic formula that relates one set of coefficients to the other. 
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